Two distinct ballistic processes in graphene at Dirac point 
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The dynamical approach is applied to ballistic transport in mesoscopic graphene samples of length 
L and contact potential U . At times shorter than both relevant time scales, the flight time tL = L/vg 
(vg - Fermi velocity) and tu = h/U, the major effect of the electric field is to create electron - hole 
pairs, i.e. causing interband transitions. In linear response this leads (for width W » L) to 
conductivity (T2 = 7r/2 e^//i. On the other hand, at times lager than the two scales the mechanism 
and value are different. It is shown that the conductivity approaches its intraband value, equal 
to the one obtained within the Landauer-Biitticker approach resulting from evanescent waves. It 
is equal to ai — A/-k e /h for W >> L and tu « tL- The interband transitions, within linear 
response, are unimportant in this limit. Between these extremes there is a crossover behaviour 
dependent on the ratio between the two time scales tL/tu- At strong electric fields (beyond linear 
reponse) the interband process dominates. The electron - hole mechanism is universal, namely does 
not depend on geometry (aspect ratio, topology of boundary conditions, properties of leads), while 
the evanescent modes mechanism depends on all of them. On basis of the results we determine, that 
while in absorption measurements and in DC transport in suspended graphene 02 was measured, 
(Ti would appear in experiments on small ballistic graphene flakes on substrate. 

PACS numbers: 72.80.Vp 73.23.Ad 81.05.uc 05.60.Gg 
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I. INTRODUCTION 

Electronic mobility in graphene, especially one sus- 
pended on leads, is extremely large [l|, so that a graphene 
sheet is one of the purest electronic systems. The relax- 
ation time of charge carriers due to scattering off im- 
purities, phonons, ripplons, etc., in suspended graphene 
samples of submicron length is so large that the trans- 
port is ballistic [2, l3| . The ballistic flight time in these 
samples can be estimated as 



tL - L/v 



a ' 



(1) 



where Vg ~ lO^m/s is the graphene velocity characteriz- 
ing the massless " ultrarelativistic" spectrum of graphene 
near Dirac points, Ek = Vg |k|, and L is the length of the 
sample that can exceed several iim [^. The extraordi- 
nary physics appears right at the Dirac point at which 
the density of states vanishes. In particular, at this point 
graphene exhibits a quasi - Ohmic behaviour, J = ctE, 
even in the purely ballistic regime. 

Determination of the value of the minimal DC conduc- 
tivity at Dirac point in the limit of zero temperature has 
undergone a period of experimental and theoretical un- 
certainty. Early on its value, a ~ 4e^//i, was measured in 
graphene on substrate [5| , yet consequently it was shown 
in experiments on suspended samples [3| that the zero 
temperature limit was not achieved, and that, in fact, 
these early samples had too many charged " puddles" , 
so that they represented an average around the neutral- 
ity point. The value in early suspended samples [2| was 
half of that and most recently settled around third of 
this value in best samples at 2K temperature [4|. Theo- 



retically several different values for the DC conductivity 
appeared. The value 



4e^ 
CTi = - — 



(2) 



had been considered as the "standard" one for several 
years [y, |7[ and appeared as a zero disorder limit of the 
self-consistent harmonic approximation [8|. It was de- 
rived for the infinite sample and this implies the assump- 
tion that the potential difference U at the contacts be- 
tween the metallic leads and the graphene flake is unim- 
portant (this is currently under intensive experimental 
[91 and theoretical [10l4l3j investigation). An alterna- 
tive and independent approach to ballistic transport in 
mesoscopic graphene samples of finite length L I l4| with 
a large contact barrier U was pioneered in [l5| follow- 
ing ideas in [l6|. They applied the Landauer - Biittiker 
formula for conductance derived for transport in (quasi) 
one-dimensional channels. 



The value 



f^2 = 771- 



77 e 
2~h 



(3) 



was obtained in the dynamical approach to an infinite 
sample [l7| and is equal to the AC value calculated un- 
der the condition oj >> T/h at finite temperatures |18l - 
l2g |. while other values like (73 = 4^ also appeared in 
the literature [2l| . The dynamical approach to transport 
was applied to the tight binding model of graphene [17| 
to resolve this " regularization ambiguity" . The ballistic 
evolution of the current density in time after a sudden 
or gradual switching on of the electric field E was eval- 
uated and approaches the large times limit a2E. The 



physical nature of the quasi - Ohmic "resistivity" with- 
out either charge carriers or dissipation in infinite sam- 
ples (namely neglecting leads or " reservoir" ) is as follows 
[171 I22I . The electric field creates electron - hole exci- 
tations in the vicinity of the Dirac points similar to the 
Landau - Zener tunneling effect in narrow gap semicon- 
ductors |23[ or electroi i - p ositron pair creation in Quan- 
tum Electrodynamics [2J, |25[ . Importantly, in graphene 
the energy gap is zero, thus the pair creation is possible 
at zero temperature and arbitrary small E, even within 
linear response. Although the absolute value of the quasi- 
particle velocity Vg cannot be altered by the electric field 
due to the " ultrarclativistic" dispersion relation, the ori- 
entation of the velocity can be influenced by the applied 
field. The electric current, ev, proportional to the pro- 
jection of the velocity v onto the direction of the elec- 
tric field is increased by the field. These two sources of 
current, namely creation of moving charges by the elec- 
tric field (polarization) and their reorientation (accelera- 
tionjare responsible for the creation of a stable current 
171 I22I |26[ . The result within linear response is that the 
current settles very fast, on the microscopic time scale of 
t^ — fi,/7 ~ 0.24 /s (7 being the hopping energy), on the 
asymptotic value. 

A deeper analysis of the "quasi - Ohmic" graphene 
system beyond the leading order in perturbation theory 
in electric field revealed 271 that on the time scale 



eEvo 



(4) 



the linear response breaks down due to intensive Landau- 
Zener-Schwingcr's (LZS) pair creation and leads to a lin- 
ear increase with timej28|. At times larger than tni the 
result is consistent with the WKB approximation [26|, |29[ 
This is in contrast to dissipative systems, in which the lin- 
ear response limit can be taken directly at infinite time. 
This perhaps is the origin of the " regularization" ambi- 
guities in graphene, since large time and small field limits 
are different. Recently the WKB approximation to the 
interband or the LZS transition was extended to the fi- 
nite samples [30[. In this case the ballistic evolution is 
"truncated" at ballistic time i^, Eq.([T]). It is (sometimes 
implicitly) assumed within both the dynamical approach 
and the Kubo approach that there is no significant con- 
tact barrier U between the leads and the graphene flake. 
In contrast, the Landauer - Biittiker (LB) approach 
hinges on the description of the leads in terms of a po- 
tential barrier of a certain non-zero barrier height U (r) 
jld |. The first quantized Weyl equation therefore was 
considered in ref. [15| in full analogy to the one-particle 
Schrodinger equation to define the T - matrix and it 
was found, surprisingly, that quasi - Ohmic behaviour 
emerges for large aspect ratios W/L — > 00 and more- 
over the conductivity is precisely ui for rectangular sam- 
ples. The approach was extended to barriers of various 



shapes [3l| , boundary conditions [ll| , nonrectangularge- 
ometries like Corbino disks using conformal mapping j32| . 
It was established that the conductance is quite sensitive 
to the topology of the sample, but rather insensitive to 
the potential shape. This approach cannot be extended 
to the AC fields. This is an additional motivation to 
extend the dynamical approach to include the effects of 
the barrier potential. The barrier potential provides an 
additional time scale 



tu = h/u. 



(5) 



It looks like the physical picture behind the LB ap- 
proach is almost "orthogonal" to the one of the ultra- 
rclativistic pair creation mentioned above. In this note 
we rigorously apply the dynamical approach to study 
transport in mesoscopic samples. We demonstrate that 
the physics behind the two values of the DC conduc- 
tivity is quite different despite the fact that numerically 
CT2 = 1.57eV/i is just 24% higher than tri = l.27e^/h for 
the stripe geometry (this is quite accidental due to the 
nonuniversality of the latter value [32| ) . A recent numer- 
ical simulation of the dynamics of a graphene ribbon [33| 
demonstrates that the value of conductance correspond- 
ing to 0-2 on a short time scale crosses over at II (via 
a series of strong oscillations) into one corresponding to 
the LB conductance in the ribbon (analogous to cti) at 
asymptotically long times. These two physical processes 
governing the ballistic transport are quite distinct. One 
is fast and homogeneous: the interband channel (valence 
and conduction "cones" of graphene near Dirac point), 
namely the electron - hole creation sometimes referred 
to (usually beyond linear response) as Landau - Zener 
tunneling, or, in particle physics, the Schwinger's pair 



creation [3J|. It is unique to graphene (and some other 
systems with similar band structure like topological insu- 
lators) and has certain surprising features. For example, 
this channel of conduction " dries out" or is depleted for 
any finite sample. The second mechanism, the intraband 
transition, despite constituting a peculiar " relativistic" 
kind of electron acceleration, is much more common. It 
is important for transport only for a sufficiently large 
contact potential between the leads and the graphene 
sample and unlike the interband channel, is a long time 
phenomenon. 

We use the dynamical approach to determine what pro- 
cess is dominant for the evolution of the I-V curve of a 
finite graphene sample directly at the neutrality point, 
with the contact barrier taken into account. We explore 
this evolution on different time scales t, which can be as- 
sociated with the frequency a; ~ 1/i for a periodic pulse 
or t he p ulse duration in the relaxation type experiments 
35I . |36[ . The physics depends essentially on the relation 
of a time scale t with respect to the three physical time 
scales tL,tu and t„/ defined above in Eqs. (|l|5|4p . It is 
demonstrated analytically that for a finite barrier poten- 
tial and finite length the infinite time limit coincides in 



linear response with a generalization of the LB calcula- 
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tion m 

We start in Section II from the definition of a model 
neglecting effects of the contact potential barrier, namely 
for tjj >> ii , tni ■ The dynamical approach to ballistic 
transport in the infinite sample is briefly outlined (this 
way a2 is obtained). The approach is generalized to the 
case of a finite sample and the decay of the interband 
channel is demonstrated. In Section III a phenomenolog- 
ical model of the graphene - lead coupling is specified and 
certain stationary properties, like the quasi one- particle 
T - matrix, are constructed. The LB results are slightly 
generalized to the case of an arbitrary barrier height U 
equal to the chemical potential of the leads. It will be 
interpreted as a long time limit of the intraband contribu- 
tion within the dynamical approach to the finite sample 
formulated in Section IV. The evolution of the current in 
graphene at Dirac point for an arbitrary potential barrier 
is given here (within linear response) as an integral. Re- 
sults of the numerical evaluation of the integral together 
with the analytically obtained short- and the long time 
asymptotes arc given in Sections V and VI for the intra- 
band and the interband contributions, respectively. The 
results are summarized and discussed in Section VI. 



II. SMALL CONTACT BARRIER: DYNAMICAL 

APPROACH TO THE INTERBAND 

TRANSITION 

A. The Hamiltonian for the infinite sample 

The electrons in a constant and homogeneous electric 
field near the "right helicity" Dirac point |6| are approxi- 
mately described by the Weyl Hamiltonian, 

(6) 



H = K + V; 

V ^ e f <fr i^lvi;^; V 



-ihvgcr ■ V; 



hvn 



A + $; 



(7) 
(8) 



where K and V are first quantized operators and the 
annihilation operator ip^ is a two component a = 1 , 2 
(pseudo) spinor. A and $ are the vector and the scalar 
potentials describing the electric field which is switched 
on at t = 0, oriented along the y axis and, importantly, 
is coordinate independent. We employ units in which 
h = Vg = 1. In momentum basis, -0^. = .^ ^ e''''''0j., 

k 

where D is an infrared cutoff and W is the width that 
also will be treated as large. When the constant electric 
field is written in a gauge that respects the translational 
symmetry, $ = 0, A = (0, —cEt), the different momenta 
decouple: 



H 



v^ -,t 



k 



^i 



A Vi 



(9) 



The model is a rather crude idealization of the experi- 
mental situation in several respects. The sample is con- 
sidered "infinite" and absolutely homogeneous. This al- 
lows a convenient use of the axial gauge invariant under 
translations. Finite width (perpendicular to the electric 
field ) generally creates no complications and the whole 
discussion can be repeated for finite W and variety of 
boundary conditions. On the other hand, finite length L 
breaks the translational invariance and simplicity is lost; 
this is one of the subjects of the present paper. The sec- 
ond idealization pertains the description of leads. It is as- 
sumed that the leads are absolutely unintrusive, namely, 
one can imagine the Corbino disc geometry or a lead with 
no contact potential difference. The more general case, 
with a potential barrier, will be treated below. 

The spectrum before the electric field is switched on is 
divided into positive and negative energy parts describing 
the valence and conduction band: 



{a- ■ k) Uk 



Uk 



-ekWk; (cr • k) Uk = ekWk; (10) 



where Zk = {k^ + iky) / e\^ is a phase and ek = |k|. 
Since all the momenta are independent in the Hamilto- 
nian Eq.® due to use of the gauge in which the electric 
field is represented via a homogeneous vector potential, 
a second quantized state is uniquely characterized by the 
"first quantized" amplitude. 



V'k(t) 



^kW 
^kW 



(12) 



which is a "spinor" in the sublattice space. It obeys the 
matrix Schrodinger equation in sublattice space: idtip^ — 
<T ■ (k-f |A) f/^k- The initial condition corresponding to 
a second quantized state at zero temperature in which 
all the negative energy states are occupied and all the 
positive energy states are empty is ip^^ (i = 0) = Wk. 

The evolution of the current density, J = —Aeipj.crij}^ , 
of a state in terms of this amplitude is 



jy{t)^~Ae Y, ^Ut)ay^^{t). (13) 

The factor 4 is due to spin and valley degeneracies of the 
Weyl fermions. To leading order in the DC electric field 
for /i = one obtains 17| a = a2, EqJSl). The same 
result is obtained for any frequency 19|, |20| . It was cal- 
culated analytically in the tight bindin g m odel for arbi- 
trary E (beyond linear response) in ref. 28| . Corrections 
to both DC and AC conductivity were computed in [27| 
and reveal that the linear response breaks down at i„; as 
mentioned above and is perhaps a source of the "regu- 
larization ambiguity" in linear response. The dynamical 
approach provides a simple interpretation for the nature 



of the excitations: the copious creation of electron - hole 
pairs or interband transitions. Beyond linear response 
these transitions can be treated either within the Landau 
- Zener (or WKB) approximation [25|, |26|, |30| or exactly 
by using Schwinger's method [27 1. 

The simple method of calculation used in the above 
works hinges on the translational invariance of both the 
sample and the electric field. However, as long as the elec- 
tric field is treated within linear response only, one may 
also consider an inhomogeneous bias. It turns out that in 
graphcnc at Dirac point, if there is no potential barrier 
like the one assumed within the LB approach (and dis- 
cussed in detail below), the current density decays with 
time even as the local electric field stays constant. 



(and A = 0), see the dashed line in Fig.l. This assump- 
tion holds even in high current experimental situations 
like the one described in ref. [30|. The current, to leading 
order in perturbation V, Eq.^ for Dirac point Ugate = 
is: 



lyit) 



-AW 



(15) 
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B. Linear response: decay of the current in a finite 
sample with a finite range of the electric field 

To model the bias voltage we assume that electric field 
is homogeneous in the segment —L/2 < y < L/2, and 
therefore can be described by a scalar potential. 



^(y) 



1 for y < -L/2 

-2y/L for - L/2 <y < L/2 

-1 for L/2 < y 



(14) 



As explained in detail in Sec. IV of ref. l27|, the current 
within the Weyl model has an ultraviolet divergence that 
should be removed in a chirally invariant manner. Since 
the present case is not different in this respect from that 
of the infinite range electric field, the details are omitted. 
After some algebra, the conductance (for large W so that 
continuum momentum can be used) takes the form 



Git) 



^^ .L- 



sin [{p - I) L/2] [Iskp - pski) sin [t {skp + £ki) 



oo J pd—0 



{p - /) Ekiekp 



Skp + £ki 



(16) 



r 



where k — l^ = Px,l — ly,P ~ Py In terms of dimen- 
sionless variables e == (skp + eki) L, S ^ L {skp - Ski) /e, 
A = L {p — I) /e, one obtains the (effective) scaled con- 
ductivity <T (t) = G [t) yy as function of time in units of 
the flight time t — 1/1^: 



a{t) 



4e^ 



°° sin (te) 



4=0 JA=a 



1 

A3 



A2 



eA 

(17) 
This function is given as the red line Fig. 2 and has the 
following behaviour. Before t = 1/2, a (t) = e^/4; 
therefore in physical units one recovers the "dynamical" 

2 2 

value CT2 = f 7^ = I^.This is just the result of pseudo- 
relativistic invariance (maximal velocity Vg) of the Weyl 
model. The effect of the finite extent of the electric field 
has no time to propagate to the center of the sample 
where the current is defined. Then the current drops fast 
and settles at t^ into a power decrease 



a t = 



e'L 

47rt 



7rt 



(18) 



In the whole range t > \/2 there is an excellent fit for 
this function: 



oil) 



1 



0-2 



7r< - 7r/2 + 1 - ^/t/2~\/l 



(19) 



Until now the linear response approximation was used. 
Hence, for a finite range of the electric field (finite dis- 
tance between the electrodes) a stationary flow state is 
only possible beyond linear response. The LZS tunneling 
over the band gap is generally a nonperturbative phe- 
nomenon. The linear response is only useful in a limited 
time, t <_tni, due to the ultrarelativistic spectrum of 
graphene [28|. The interband processes beyond i„; are 
discussed next. 



C. The LZS conductance beyond linear response 

Since there are two characteristic times beyond linear 
response, II and i„;, three possibilities exist: II >> tni, 
t^ « tni or II ^ tni- In most transport experiments the 
ratio t^/tni — L^JeE/hvg is smaller than 1. The time 



scale on which nonlinear effects become dominant is not 
always very large; for example, for experiments not neces- 
sarily dedicated to large current transport measurements 
37 1 in which E = W^V/m and L = 0.3/i77i, nonlinearity 
sets in at i„i — 0.3ps, which is of order of the ballistic 
time. Moreover, graphene flakes under large fields of or- 
der 2 • lO^V/m have been recently studied in specially 
designed high current density experiments [30|. For a 
sample length L = 2fxm this results in a very large ratio 
tb/tni = 100. Therefore it is of importance to calculate 
the conductance for arbitrary tL/tni- We start from large 
fields for which the LZS tunneling is the most effective. 



th » tnt 



Interestingly this formula is very accurate already at t 



tni, see ref. 28 



The physics of pair creation is highly non-perturbative 
and non-linear in nature and therefore, instead of the 
linear response, Schwinger had to use functional methods 
to get an exact asymptotic formula. The rate can be 
intuitively understood using the much simpler instanton 
approach originally proposed in the context of p article 
physics [25| (extended later to low dimensions 12J]), that 
is known in condensed matter phy sics as the Landau - 
Zener tunneling probability 2Q, [SOJ, l38| . 



The density in the infinite sample, calculated using the 
simple Landau - Zener creation rate expression for one of 
the flavours, is |25l.l26|: 



Analytic and numerical solutions of the tight binding 
model[28|, as well as of the Weyl model describing the 
physics near the Dirac point demonstrated [26|,|27| that at 
tni the creation of electron - hole pairs becomes dominant 
and is well described by an adaptation of the well - known 
Schwinger electron - positron pair creation rate -jiN oc 

{eE) ' . The difference with the original derivation [3J| 
in the context of particle physics is that the fermions are 
2-1-1 dimensional and "massless", thus magnifying the 
effect. The polarization current is J (t) = —2evgN{t) 
and therefore Schwinger's creation rate at asymptotically 
long times leads to a linear increase with time: 



a (t) ^ a2 {eEfh. 



(20) 



Ni, (t) = e{ky)e {eEt - ky) exp {-irkl/eE) , 



(21) 



where Q are the Heaviside functions. One considers the 
"tunneling" from the conduction band to the valence 
band at fixed k^. The gap is given by 2|fc2;|. Conse- 
quently the number of pairs is 



N{t) 






t, 



(22) 



(27r)^ 
corresponding to the Schwinger's pair creation rate 

■^N = -^ (eE) ' . There is a simple relation within the 
Weyl model between the rate and the current density, as 
was shown recently [26| (see also [2J|): 



J 



2e 



■7r2 , 



eEt 



kl + {eEt -kyY 



=iVk (t) 



kl + {eEt - ky) 



eE 



dt 



Ni,{t) 



(23) 



Substituting Ea. ([2T|) . one obtains: 



^<"^sr.(v^ 



+ {eEtf + k^]exp[ ^ 



7d^ 

'Te 



(24) 
At large f, one arrives at conductivity. Eg. ([20]) . 

Adaptation of the instanton approach to finite length 
sample is quite cumbersome, however the long time limit 
is simple, as was shown in ref. |30J . The range of inte- 
gration over momenta is determined by the semiclassical 
condition for tunneling: 



2e 

J{t^ (X)) = ^T 



eEL/2 fk 



fex=0 



=0 ^k 
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(25) 



where fc™^'' 

ity, 



{eEL) — k^.. The result for conductiv- 
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(26) 
-1 , 



is presented as a green line in Fig. 2. For ti » tni it is 
proportional to the ratio t^/tni '■ 



e'L 



tL 



cr (i ^ oo) = — 5-ve£/ = -^02- — ■ 

TT^ TT^ tni 



(27) 



This is larger than ui for thjtni > 7r^/4 « 2.5. For these 
fields the slow decrease due to finite extent of electric 
field, L, does not materialize and increases monotonically 
as function of time t (physically representing the pulse 



U(y) 
V(y) 



VJ2 




^ -V/2 
V o_ 



FIG. 1: Potential barrier U(y) (green line) describes contacts 
of the leads with the sample, while the bias electric potential 
V(y) (blue dashed line) describes the applied constant electric 
field (from -L/2 to L/2). 
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t (in units of ti) 



FIG. 2: Evolution of the scaled conductivity, Eas. (|3ll9l26l27p 
at Dirac point for a finite range electric field (and no contact 
barrier) . Electric field strength determines the time scale f „; , 
Eq.(|4]) at which nonlinear effects set in. 



period or pulse duration), 
situation is different. 



For small ratios tL/tni the 



In this case of small electric fields one first encounters 
at II a powerwise drop in current density from the short 
time value (T2 as in Eq. p^ before the nonlinear effects 
take over. For large times t > tni one can use the same 
semiclassical LZ method to determine how the power- 
wise decrease is halted. The asymptotic value given by 
Ea. (l26l) is this case simplifies into: 



2e^ 



CT (i ^ oo) = ^^ (eE) L^ = -^CF2 — 

T^ T^ \tnl 



tL 



(28) 



This intercept with the slow decreasing current, Eg. psp . 

see Fig. 2, occurs at very large time t = f "f^- When t^ 
is of the same order as t„/, the effective conductivity is 
approximately (72, see values close to tL/t„i = 7r^/4. 

Till now we neglected the influence of a potential bar- 
rier between the leads and graphene on the ballistic trans- 
port. 



III. CONTACT BARRIER: STATIONARY 

PROPERTIES AND THE FIRST QUANTIZED 

LANDAUER - BUTTIKER APPROACH 

A. Phenomenological description of contacts. 
Symmetry of the Hamiltonian. 



U{y) 



for y > L/2 or y < -L/2 
[/ > for -L/2<y< L/2 



(29) 



The derivation of the model from a microscopic Hamil- 
tonian is discussed in several works, see for example [lO| 
and it was found to describe a typical transport experi- 
ment quite well. 

The second quantized Hamiltonian is H ~ K + V , with 
the perturbation (bias) given in Eq.([8]) and the modified 
"large" part 



K 



tprHlQtpr 



(30) 



The "first quantized" operator Hiq now contains the bar- 
rier potential U (r): 



H 



iQ 



V + C/(y) 



(31) 



The barrier breaks the translational symmetry, however, 
for the simple form of the symmetric barrier we have 
adopted, Ea. (|29p . the operator Hiq is invariant under 
reflection, P : y -^ — y, supplemented by the spinor rota- 
tion. 



S^Pa, 



[HiQ,S]=0. 



(32) 



The bias potential, Eg. p^ . is also chosen to be antisym- 
metric which simplifies the considerations. 

The presence of a barrier renders the problem analo- 
gous to that in mesoscopic physics [23| and suggests that 
the T - matrix approach to transport is useful in this 
case. 



One models the effect of coupling to leads by a finite 
(and sometimes very large 15|) potential energy barrier 
(that should in principle be found self consistently jlO| ) . 
The simplest model is the square barrier, see Fig.l, 



B. T - matrix 



The LB approach utilizes the notion of a transmission 
coefficient through the channel n, T„ = |t„ (^)| , where 



tn is its amplitude. The conductance is 



g(^) = tE^- 



(33) 



where the summation is over all the open channels. 
Therefore one should solve the " classical" Weyl equation 
with a barrier 



HlQ^ 



W + U{y)]ijj = eip. 



(34) 



For simplicity we consider only periodic boundary con- 
ditions in the direction perpendicular to the field with 
"perimeter" W, although various more realistic bound- 
ary conditions were discussed in ref. [I3] and numerous 
works since. Due to translational symmetry in the di- 
rection perpendicular to the barrier, we consider a fixed 
value of the momentum p^. = fc = -^nx- Despite the lack 
of translational symmetry in the field direction y due to 
the barrier, one can still use the momentum py = p as 
a good quantum number for scattering states. Another 
" number" is the sign of energy which determines the wave 
function in the leads, namely distinguishes between the 
u and the v spinors given in Eq. (fTTj) . The reflection sym- 
metry defined in the previous Subsection converts left 
movers into right movers 



^'^k,p^ — ^p'^k.—p^ 5 ^'^k,p^ — ^p'^k,—p^ 5 

(35) 
where we suppressed the index k in Zkp = {k + ip) /ck- 

The "out of barrier" equation is just the free Weyl 
equation with negative and positive energy solutions, 
^ = Ukpe'^''''+Py^ (hole) and i' = Wfcpe*('''^+™) (electron) 
discussed in the previous Section. It should be matched 
with the " in barrier" solution. Several distinct kinematic 
possibilities exist. We survey them from high to low, see 
Fig. 3. 

1. Energies above the barrier, e > U. Both inside and 
outside one has electron v- states with different momenta. 
Outside the barrier p = y/e^ — fc^, while inside the barrier 
momentum in the field direction is 



One has a wave (real q) inside for 



P>P^= v/(C/ + 2|fc|)[/. 



(36) 



(37) 



The lower bound, pi (fc), is the green line in Fig. 3. There 
is an evanescent (imaginary momentum q) particle state 
inside for 



\/C/2 - fc2 =pu{k) <p< P2. 



(38) 



Crossing the red line in Fig. 3, pu (fc), one encounters 
states below the barrier. 



2. Positive energy states below the barrier, < e < U. 
One has the v spinor (electron) outside the barrier, while 
the u spinor (hole) inside. For momenta p in the range 



^iU-2\k\)U = pi<p<pu, 



(39) 



the states are evanescent hole states. At yet lower ener- 
gies, p < pi, one has a propagating state, but this time 
a hole. This relativistic feature is the cause of the Klein 
paradox. 

3. Negative energy states, e < 0. Outside the barrier 
now one has e = ~ \/pP^+\? . This is another purely " rel- 
ativistic" possibility in which one has holes both outside 
and hence inside the barrier. 

The Schrodinger equation above the barrier e > C/ is 
solved by the scattering states for right movers, p > 0, 

*.M-^, (40) 

Vkpe^'^y + rupVk.^pe-'^'Py for y < -L/2 
X { Akpykqe'^y + BkpVk,-qe-"'y for -L/2<y< L/2 . 



tkpVkpe'Py for L/2 < y 



Matching conditions. 



tpVpe'P"^^^ = ApVge^^'^/^ + BpV-ge-'^'^^^, 

(41) 

(where the fixed index k is suppressed) determine the T- 
matrix and are easily solved. The electron v states have 
to be replaced with u states in the case of a hole, so in the 
second energy region in the barrier part v -^ u, while in 
the third energy region in all parts v -^ u. For example, 
for evanescent modes below the barrier one obtains 



tkv 



^ — iqL 



(1 - ZpZq) 



piqL 



(Zp 



(42) 



C. LB conductance 

According to the LB formula one has to sum up 
the transmission coefficients Tkp = \tkp\ over all the 
states, namely over the regions 1-3 that obey a constraint 
p (k) ~ \J [i^ — k^. The chemical potential fi ~ U + Ugate 
of graphene is counted here from the bottom of the bar- 
rier and let us consider only positive Ugate- In this case 
region 3, in which energy is negative, and region 2, where 
energy is positive but not large enough therefore will not 
contribute and we are left with region 1 of the previous 
Subsection: 



9-77 /-^ 



2tt 



'^k,p{k)\ 



(43) 
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FIG. 3: Kinematics of scattering states of the first quantized 
Weyl equation with potential barrier. The red line corre- 
sponds to p,l = \/U^ — fc2, the green to p, / = y/U {U + 2k), 
and the blue to p, / = yJU {U — 2k). The lines separate kine- 
matical regions for the intraband transitions. Here U = 1. 



the electron - hole "bulk" dynamics, to the long ballistic 
time, the barrier reflection dominated dynamics. 



D. Electron tunneling into graphene from leads for 

U = Ugate 

The main physical effect of leads is that they " contam- 
inate" the graphene flake (assumed to be at Dirac point). 
Electrons from the lead metal tunnel into the flake cre- 
ating charge "puddles" on both sides of the barrier. As 
a result the intraband channel for electric transport is 
greatly enhanced for small flakes and might dominate 
over the interband channel, especially at small fields. The 
quantitative consideration of this contamination is more 
transparent in the basis of symmetric and antisymmetric 
eigenfunctions than it would be in the scattering states 
basis (which was convenient for the LB approach above). 

The symmetry image of the scattering state (j)^, (r) 

defined in Eq.(j40|) is 4>k,-p{'i') = S(j)^.p{r), where the 
symmetry operator S was defined in Ea. (|32p . Above 
the barrier (region 1), e > U , one can write symmet- 
ric (s)/antisymmctric (a) functions as (with the x de- 
pendence of the wave functions e**^^ implied and index k 
suppressed) : 



Let us first consider, following ref.[l5[, only evanescent 
states contributions under the barrier. Although finite 
sample width with various boundary conditions can be 
easily considered, we take a limit of infinite aspect ratio 
W/ L — >■ oo and generally replace summation over k = 
kx in Ea. (|43)) by an integral. In this case one can form 
two dimcnsionless combinations: the barrier "strength" 
UL = fl and UgateL. It turns out (not shown in the 
present paper) that the conductance has a smooth limit 
UgateL — > that, of course, involves evanescent states 
only. One therefore can study directly the case of Ugate = 
0. Substituting the transmission coefficient of Ea.(|42|). 
one obtains the following limiting value of conductivity: 



W 



u^-e 






2e^VL 



1-t 



71-2 j-^^ ^^^j^2 /^^N 



r 



(44) 



where the dimcnsionless momentum k = k/U is used. 
The result for various J7 appear in Fig. 4 as the long time 
limit. 

For large J7 the integral is dominated by small k and 
one gets the "mesoscopic" value of conductivity, 
»i 



CLB {^) 



2e'n 



2e' 



k=0 



cosh~2 (fcO) =: :l^ = cTi. (45) 



One therefore can apply the dynamical approach to try 
to understand the crossover from the short ballistic time. 



1 



t.ri,)^-^^(l±S),,(y)^^^_ 



1 



(46) 



{e'Py + a^Dl^'^e-'Py) Vp for y < -L/2 

± {e"^y + <Jxe-"^y) E^'^^Vg for ^L/2<y < L/2 , 

± (D'pe^Py + a^e-^Py) Vp for L/2 < y 

where q is the momentum in the barrier range given in 
Eq.dSSl), and 

'ipL ^p^q ~ 1 =1= (Zp ~ Zq)e 



or = T^'p 



Zp- Zq± {ZpZq ~ 1) e"'?-^ ' 



-l{q+p)L/2 



E^'- = Ap±BpZ_q 



{4 - 1) 



Zq ± {ZpZq - 1) e-*?i ■ 



Below the barrier for the positive energy states (region 
2) one similarly has the same form 



{y) 



1 



y/2DW 



(48) 



{e^Py + a^D'^'^e-'Py) Vp for y < -L/2 
X < {±e'''y + a^e-'iy) E'^^°-Uq for -L/2<y< L/2 
± [D'^-^e^Py + a^e-^Py) Vp for L/2 < y 

with different coefficients 

rt's,a _ -JLp 1 + ^p^g "F [Zp + Zq) e 



+ ZqT{zpZq-l)e-iLi'' 



'iL{p+q)/2 



Ep 



.4 - 1) 



Zq + ZgT {zpZg + 1) e-ii^i ' ■ 



The negative energy states (region 3) are obtained from 
those in region 1 by a replacement u — >■ u. 



vTiy) 



1 



^ ^'" V2DW 

e'Py + a^D^p^^e-'Py) Up for y < -L/2 



± {e"iy + a^e^'iy) E^'^Uq for -L/2<y< L/2 . 
± {D'pe'Py + a^e-'Py) Up for L/2 < y 

We infer that physically the leads induce charge "pud- 
dles" in graphene at Dirac point. These electrons can 
now be accelerated (in an ultrarelativistic fashion by re- 
orientation) and compete with the electron - hole channel 
described above. 

Now we calculate the evolution of the current density 
within the linear response approximation. 



IV. CURRENT EVOLUTION IN GRAPHENE 
WITH BARRIER 



where 



i;-+ST 



.■5t 



viy)^tHy)rkp(y) 



(55) 



(50) In the middle of the sample the current operator is 



4(y-o)= jy{x,o) 



(56) 



where, for example, 

ji'^^M = wskk'jt;^; 3tr/' - -4e</>rj (0) <Ty^ii (0) . 

(57) 
The factor 4 appears due to spin degeneracy and two 
Dirac points (valley degeneracy). Now one can write the 
first order contributions to the electric current induced 
by the perturbation. 



A. The bias and the electric current in basis of the 
barrier eigenstates 



B. General expressions for the interband and the 
intraband contributions to current in linear response 



The dynamics of the manybody system with the con- 
tact potential barrier after switching on the electric field 
is determined in the Heisenbcrg picture by the equation 



idt^ (t) = [K + V)i; (t) , 



(51) 



In linear response one obtains at gate potential Ugate 
(assumed positive) two contributions with completely dif- 
ferent physical interpretations. In the first term the sum- 
mation is over electron states above fj, ~ U + Ugate and 
electron states below the Fermi level. 



where the K is defined in Eq. ((SO]) and V in Eqs. ((8l[T4l). 
The unperturbed part of the Hamiltonian, A", is diag- 
onalizcd by choosing a basis linked to the eigenfunctions 
of the first quantized operator Hiq, Eq.([34|): 



lee (t) 

-^E E E 

k p-ekp>t^l:ekl<f^ 



I _ g-'»(efcp-EA:()t 

£kp — £kl 



+ST-++TS 



'^kdp Jk 



(58) 



cc. 



''kp 



i'r 



SO that K 



A^t 



= 4>Vp(r)^r-Mp 



s f^\ is 



(52) 



E [^kp (r) alp + ^tp (r) hy,p 

kpS 



■J. 



bibx 



where A combines 



Ea £fcp [a'xa\ 

k,p,S. The operator a^ (b) creates an electron (hole), 
and summations are over symmetric and antisymmetric 
states, S = a, s, with p > 0. The perturbation operator 
in the new basis takes the form 



v^T v++^' 



X,a 



- + ^t. 



■Z,-^a{a, + V^/h\a, + V^' a\h, + V^-blb, 



(53) 

For example, the electron - hole matrix element is given 
by 



V, 



hST 



kl-,k' p 



<''-'>V{y)^l}{y)4>l,{y) 



W5,,yV, 



VST 
kk> ^k,lp 1 

(54) 



The summation over symmetry indices T and S is un- 
derstood. The second contribution sums over electrons 
above ^ and all the hole states with arbitrary I: 



leh (t) 

k,l p:ekp>IJ- 



(59) 



-iiskp+i^kl)t 



£kp + Ski 



'^k,lp Jk,pl ^ ^^■ 



The hole's momentum has no restriction since its energy 
is always negative, while Ugate is positive (in this work we 
do not consider the case Ugate < , that was discussed 
recently in ref.|3l| in the framework of LB approach, 
where other possibilities can occur). The first contribu- 
tion is the "one- particle" type (the intraband channel), 
very much like in more common manybody electronic 
systems. The second contribution, to the contrary, is 
purely ultrarelativistic (the interband channel) and de- 
scribes the electron-hole pair creation, very much like in 
the infinitely long flake discussed in Sec. II. 
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It is useful to define an efi^ective dimensionless conduc- 
tivity. There are only two time scales in the problem: 
the inverse barrier height t(j = 1/U [h/U in physical 
units) and the flight time to cross the sample ti, ^ L 
(L/vg) with the ratio being the only "material" param- 
eter n = UL ^ t'L/tij in the linear response considera- 
tions. For this case the effective conductivity has a scal- 
ing property: 



a{U,L,t) = 



I{U,L,t) L 
Vo W 



nj) 



(60) 



The scaled time is defined, as in Sec. II, t ~ t/t^. Since 
in the dynamical approach the limit Ugate — > is smooth 
[17|, in what follows the conductivity will be calculated 
directly for Ugate — O.We concentrate on the case when 
quantizations of k^ and ky are not important (the gen- 
eralization to the discrete case is straight forward), the 
sums in Eqs. (|58|) and (|59|) can be replaced by integrals 
over dimensionless momentum variables, k — Uk, p — Up 



and pu 



a- (t) 



2 

1 — k . The intraband contribution becomes 



LU^ W f D 



Vf) 27r V27r 



k—Q J l—O ^ p—Pu 



^Ip Jpl 



(61) 



cc. 



Here rescaled energies, e = Skp/U, e' = eu/U in addi- 
tion to momenta k,p,l are used and the factor 2 comes 
from positive and negative k. Similarly the interband 
contribution to conductivity takes the form 



for the electron - electron contribution for T ^ S, 
ieVo 



V, 



++TS 
Ip 



X Ul + z*iZp) 



2DW 

jjiT*jjS^t{p~l)L/2 _ g-j(p-0L/2 



(63) 



p~l 



j^S^i(p+l)L/2 _ jjiT*^-i(p+l)L/2 

H^P^^!)^ ^T7 

^AE[^*EI [(1 - Zgz;) f {g-q) + {z, - z;) f{g + q)] } 

where the function / is defined as 

sin iqL/2) L cos {qL/2) 



/(?) 



(64) 



It is important to note that the dominant contributions 
to the integrals in Ea. (|6T|) come from the {p — l)~ poles. 
The corresponding current density matrix elements are: 



— \-sa 



■'P' WD 



E;*E[''{l + z;){l + Zg); (65) 
-E"CE','' (1 - z!) (1 - z„) . 



The one-particle (electron - electron) contribution to 
the scaled conductivity. Eg. (150]) is shown in Fig. 4 as 
red curves for various values of 51 as function of time t. 
At times shorter than both t^ and tjj it rises linearly. 
For 51 > 1 the conductivity oscillates before eventually 
approaching the LB result of Sec. Ill C. In the following 
subsections the two ballistic regimes, namely, the limit 
of long times t » tjj, tL (or in scaled form 1. » 1/Sl, 1) 
and short times t « ic/i^L (or in scaled form t << 
1/51, 1) are analyzed analytically. 



^eh 



t) ^~2- 



LU^ W f D 



k=0 J 1=0 Jp=Pu 



Vn 27r V27r^ 

1 _ g-i(e+e')^^* 



(62) 



-V, 



-+ST ■+-TS 



Ip 



Jpl 



Next we consider the two contributions in turn. 



B. Short time asymptotics 

Expanding the electron - electron contribution to the 
conductivity, Eq.(|6T|). to leading order in the short time 
limit yields 



V. THE INTRABAND CONTRIBUTION 



A. The electron 



electron matrix elements in linear 
response 



Neglecting the dependencies on the cutoff and keeping 
notations used in Sec. IV A one obtains, after some al- 
gebra, the following expressions for the matrix elements 



= -it- 



-e^n 



Pu 



fc— 0"'/— <J p^Pij 



^Ip 



V,r"'j.r''+cc. (66) 



This leads to a small conductivity contribution since the 
integration range on the variable I is rather limited. The 
linear behaviour of the conductivity is shown in Fig. 4 
(red lines). For the case fl < 1 [II < tu), represented 
in Fig. 4 by 51 = 7r/16, 7r/4, the intraband contribution 
is positive and increases monotonically. However, when 
51 > 1 {tL > tu) represented in Fig. 4 by 51 = tt, 27r, it 
becomes negative. 



C. Long time asymptotics 
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Due to the oscillating functions in Eq. (|6T|) the long 
time asymptotics is due solely to the region of the three 
dimensional integral when ekp — Ski — > 0. Consequently, 
in view of the discussion of the various kinematical re- 
gions in Sec. Ill B, summarized in Fig. 3, the limit is dom- 
inated by integrating over the transitions from evanescent 
states above the barrier (region 1) to evanescent states 
below the barrier (region 2). Since the integral is dom- 
inated by the neighbourhood of the line ekp,Ski ~ U, 
one can expand around this line (the red line in Fig. 3), 
P ~ Pu + ^P' ^o that the evanescent momentum (in the 

—2 

field direction) is, to second order, k — ^Ap^. Simi- 
lar expansions are made for the momentum I. Therefore 
one can simplify the expression for the conductivity by 
replacing the integrand in Eq. (|6ip by its limiting value. 
Albeit, this limit is nontrivial. In particular, the lower 
components of the spinor in Eq. lfTTj) is just 



Zp^Z^k + ipu, Zq = --^Ap. 

2k 



(67) 



As a result the coefficients determining the eigenstates, 
Eas. (|47l49p . simplify considerably: 



rpa.s 



P 



j-ifa.s 

El' ^e 



n{k~tp„)/2_ 



1 



D 



— > -e 



-iQp,- 



1 + Ze"^ 



(68) 



Z±e 



Vlk 



Since at small k, / (k) :^ -^k, the third contribution to 
the perturbation matrix elements, Eq.(|63|). is negligible. 
The results for the matrix elements of the perturbation 
and the current are 



Vas '\;sa* 
Ip ~ ^ ^Ip 



■sa 
Jpl 



2eVU 



Pu 



DW 

8ie 



ipij — sinh (flk) p ~ I 



(69) 



pI 



DW ip^j 



sinh (J7fc) 



Substituting these matrix elements into Eg. (|6T|) . one ob- 
tains the contribution to conductance: 



a- (t) = 



16e2f7 



pI (k) 



Ti" ./fc=o cosh (Qk) - k 



(70) 



Pu r°° 



sin [(e - e') nt] 1 



Since p 
integrals 



/=0 Jp=Pu £ ~ ^ P— i 

e — e') /pjj, one gets for the last two 



Pu(k) 
1 



sin [(e - e') fit] 



k Je=l 



{e~e'Y 



.- _ ,,,s{k,^) 
Pir{k) 

where the function 

nl 



(fc,- 



sin [(e — e') r] 



(71) 



(72) 



Therefore at large times t > tjj, one finally obtains 



cr'^" (ft) 



2e^n 



l^t 



=0 cosh^ iyik) - k 



<JLBm- (73) 



This is one of the main results of the paper. The elec- 
tron - electron contribution converges to the LB result, 
Eq.(|33]) at large times. Furthermore, with large J7 the 
conductivity is dominated by infinitesimal k and conse- 
quently can be written as 



s^r!^ '•1 



e^ 1 

^Z- (74) 



71-3 2 77:^0 COSh^ {kn) 2-K TT 

Thus, for four fiavours, we reproduce cti, Eq.([2]), starting 
from the dynamical approach . 



VI. THE INTERBAND CONTRIBUTION 

A. The electron - hole matrix elements in linear 
response 

Similarly to the above, the matrix elements for the elec- 
tron - hole Landau-Zener-Schwinger contribution read 

(with T 7^5-): 



V, 



-+TS 



Ip 



ieV 
2DW 



(75) 



jjS^i(p+l)L/2 ^ jj*S^-i{p+l)L/2 
jj*S-qS^i(p-1)L/2 ^ ^-i(p-l)L/2 



[1 ~ ZpZi 



'I ^P' 



p~l 



+AEf*El [(1 - z,z-^) f [g -q) + (z, - z'^) f (g + q)] } 

Unlike the intraband contribution of the previous Section 
the p — I pole is cancelled by a numerator (both p and 
/ are positive) and therefore the contribution will not be 
dominated by the evanescent states. The correspondent 
current density matrix elements are 
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tH — sa 



„-H — as 
Jpl 



^E;*Et{l + z~'){l + z,) (76) 



1 P 



1 



')■ 



The expression for conductivity a'^'^ {U, L, t) is UV 
divergent like the conductivity of the infinite sample 
(T (C/ =: 0, L, t) = (To (f) biased in the region of length 
L that was studied in Sec. II. B (and which is solely due 
to the electron - hole pairs). The difference Act^'' (<) = 
a'^^ {U, L, t) — (To (i) however is finite. Using Eq.([62|) for 
(T^'', the difference can be written as an integral over a 
limited domain of momenta in the field direction: 



0.25 
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n = ;r/16 





=~«=^ 





0.45 
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2 4 6 8 10 
t (in units of ti) 
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A(T 



eh 



t) 



nu{WDy 

47r3Vb 

Pu I _ g-i(e+£')nt 



/=0 Jk=0 Jp=0 



(77) 



v,r""j,r""+-- 



The results, for Ugate = 0, are given in Fig. 4 for several 
n as blue lines and were discussed in Sec. II. At small 
times it starts with the ultrarelativistic value cr2 = ^ and 
at relatively large times (% » t > t^^) it decays as -i=. 
This short time value (T2 does not change when U > 
provided the time is smaller than iL/2. This follows from 
the fact that in relativistic graphene information about 
barrier cannot arrive at the center of the sample before 
that time. The long time asymptotics in the presence of 
the barrier is considered next. 



B. Long time asymptotics 

Due to oscillations the long time behaviour of the 
electron - hole contribution is dominated by the region 
£kp + Eki ^ in the integrations in Eq. ([77)) . Unlike for 
the intraband contribution in the last Section, in this 
case all three momenta k,p,l in Ea. (j77p are small. In 
this limit q,g — >■ 1, and hence the phases Zq,Zg — >■ i 
and the function f {g + q) -^ — ^ [ficos (fi) — sin(ri)] 
and f {g — q) — >■ 0. Therefore it is simple to calculate 
Acr'^'' for special values of il. This is done in Appendix. 
Numerical results show that Acr^'* (f^ = 37r/2 + 27r A^) , 
Acr^'' (7r/2 + 27rA^) oc 1/f, while Aa^'' (tt + 2ttN) = 
—Aa^'^ {2irN) « l/47rt. One can fit the long time asymp- 
totics as Aa^'^~ cos (r2) /Ant. 



VII. DISCUSSION AND CONCLUSIONS 

Two different kinds of ballistic behaviour occur in un- 
doped graphene (graphene at Dirac point) at zero tem- 
perature. The first one is a very unusual "ultra - relativis- 
tic" interband physics. Electron - hole pairs are copiously 



FIG. 4: Conductivity of finite length samples for Q = 
UL/hvg = 7r/16, 7r/4, TT, 27r. The intraband contribution is 
the red line, interband contribution - blue lines. The total 
conductivity - the black lines. The time is given in units of 
tL ~ L/vg, while unit of conductivity is e^/h. 



created via Landau-Zener-Schwinger's mechanism by an 
applied electric field. It is not dependent on leads and fi- 
nite size effects of the graphene sample. On the contrary, 
the second one, the intraband physics is mostly sensitive 
to finite size effects and contacts. The interband transi- 
tion results, within linear response, in the universal bulk 
value (T2 of conductivity, while the intraband transition is 
characterized by a shape dependent linear response with 
the effective conductivity ai for large aspect ratio rect- 
angular flakes. 

Now we recapitulate under what conditions either of 
these two processes is dominant in experiments on a time 
scale 1/oj in an AC electric field E (or in a pulse of du- 
ration 1/uj) for a graphene flake of length L (in the limit 
of large width W, although the discussion can be easily 
extended to finite aspect ratios) and a contact barrier po- 
tential U. Here we classify various practically important 
ranges of sample (U^L) and experimental {u),E) param- 
eters. 



A. "Unintrusive" experiments, (7 = 

In reflectance and transmission experiments in visible 
to mid IR or even microwave range [39[ there are no 
leads, hence no potential barrier, [7 = 0. The inter- 
band (Landau - Zener - Schwinger) process is dominant 
for any practically length and electric field E. However, 
the transport can be either linear or highly nonlinear. 

(i) For l/w < tni = yJh/eEvg one has linear response, 

J = (J2E, with the interband value of conductivity (T2 = 

2 
■j^. In this case the nonlinear Schwinger 's pair creation 

regime is not yet reached. The conductivity is real and 



frequency independent (no inductive part). There are 
small frequency dependent corrections to linear response 
(the E^ correction to both real and imaginary part and 
third harmonic generation was calculated in ref . j27l | , see 
Eqs.(53), (54) therein). 

(ii) For 1/uj > tni and tni < t^ = L/vg the transport is 
still dominated by electron - hole channel, but is nonlin- 
ear. The electron - hole pairs are efficiently created due 
to the LZS mechanism with rate proportional to E^/"^. 
This results in following I - V curve, see Eg. ((27)) : 



7r2 V ft 



(78) 



(iii) For l/io > tni o,nd tni > ih the transport is still 
dominated by electron - hole channel and the LZS pro- 
cess but since the electric field is applied in the limited 
space (length L which is not large enough) and the cur- 
rent is much smaller, see Eq. ([28| : 



(79) 






eE 
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the sample. There is no effect of ballistic acceleration 
across the sample for large frequencies or short pulses. 



C. Small barrier 

For certain materials and geometries the potential bar- 
rier may be much smaller, of the order O.Oley (for exam- 
ple the Ti lead in the /3 = 7r/2 geometry has —0.04, while 
Ni in geometry j3 = tt has 0.04ey, see reference [13[ for 
details). It is feasible that the barrier even can be tuned 
to zero. It is not clear what is the contact barrier in 
the suspended graphene systems. It is reasonable to as- 
sume that for suspended samples the barrier is smaller 
due weaker coupling. Moreover, even in two probe exper- 
iments on suspended graphene the Dirac point appeared 
exactly at zero bias voltage [2|,|3| (compared to up to 40^ 
in sampled on substrate) and there are less charge " pud- 
dles" . That also signals that the leads do not create the 
puddles. In such a case a variety of crossover phenomena 
might occur when t[/ is comparable to t„/ or t^ studied 
here. 



D. Finite size and topology effects 



B. Large barrier 

In samples on substrate with metallic leads the work 
function of the graphene and the metal is typically dif- 
ferent and as a result the contact potential difference is 
of order U = 0.1 — leV , see calculations in jl3[ and ref- 
erences therein. It can be both positive and negative. In 
this case the corresponding time scale tu ~ h/U < 7fs 
and typically smaller than any of the other scales t^ ~ 
L/vg, tni = y^h/eEvg. This leads to an effective suppres- 
sion of the electron - hole channel for all the frequencies 
in the infrared range and smaller (including DC) and the 
physics is dominated by the electron - electron channel. 

(i) l/o; > tjj. The DC conductance is given by the Lan- 
dauer - Biitticker formula, generalized in Eq. (|44p , and is 
more sensitive to the properties of the leads than those of 
graphene. When graphene is "nominally" at Dirac point, 
namely, when the chemical potential of the lead is on the 
barrier, graphene is still contaminated by charges tun- 
neling into the stripe from the leads. These electrons are 
accelerated and lead to the mesoscopic type of conduc- 
tance. For a large aspect ratio the effective conductivity 
is ai = - 7- . The assumption of an " infinite" barrier was 
made early on in |15| in order to develop the mesoscopic 
approach to transport in graphene. 

(ii) l/w < tij. The high frequencies (microwave and 
above) experiments are done without leads. However if 
one had a set-up with leads it could not significantly al- 
ter the pseudo-Ohmic behaviour with a = a2 since the 
contaminated regions constitute only a small fraction of 



The calculations made in the present paper can be 

easily generalized to a finite or even small flake width 

W (assumed infinite here). In fact, it was shown 

numerically [33[, that in graphene nanoribbons modeled 

by a tight binding model one observes a crossover from 

2 
(T2 to airibbon = x ^^ ^ time roughly corresponding to 

^L- The value of (Tiribbon in this model is different from 
^y due to the different topology and was calculated in 
j40| . In Fig. 3 of [33 the time evolution of current for 
i = — lO'^ti is shown. The current rises from zero and 
settles on CT2(via oscillations) on the microscopic time 
scale of t^ = fi,/7 « 2.6 • 10~^^s. This time scale is not 
seen in the Weyl model (it appears as an ultraviolet cut- 
off). Yet similar results, including oscillations, were ob- 
tained in the solution of the tight binding model of the 
infinite sample [17|. The crossover from the initial inter- 
band to the intraband behaviour also occurs as a series 
of oscillations around the value uiribbon- The period of 
oscillations in Fig. 3 in reference [33[ is roughly i^, as it 
is also the case in our Fig. 4 for large D.. It is not easy 
to compare these results directly, since the microscopic 
model for leads cannot be readily translated into a poten- 
tial barrier U of the continuum model. However, the sim- 
ilarity with our results suggests that fi = UL/Hvg >> 1, 
since their number of oscillations is large. 

Our calculation is trivially extended to any system 
with Dirac point - like spectrum like double layer 
graphene and recently synthesized family of materials 
called " topological insulators" 



411 in which a surface ex- 



citations are similar to those in graphene with notable 
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exception of the chirality. Schwinger's mechanism is also 
expected in these materials since chirality was not in- 
volved (left and right movers contributed equally to the 
emission rate of graphene) . On the contrary, the effects 
of the sample topology probably require applications of 
sophisticated methods like conformal mapping used in 

Finally, let us remark on the role of the Klein paradox 
in ballistic transport in graphene. While the pair creating 
mechanism can be linked to the Klein paradox behaviour, 
the intraband mechanism at Dirac point does not make 
use of the states undergoing the Klein paradox. The 
Klein tunneling states described in Sec. Ill B (2) are 
characterized by energies in the range < e < U, the 
states below the blue line in Fig. 3. These states make a 
contribution to the current at finite times, however their 
contribution to current at large times vanishes. 
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APPENDIX. LONG TIME ASYMPTOTICS OF 

THE INTERBAND CONTRIBUTION FOR 

INTEGER VALUES OF fi/ (27r) 



In the long time limit of the e-h conductivity, Eq. ([77| , 
the "under the barrier' momenta have a simple limit: 



1, so that the phases Zq,Zg — )■ i and the function 

0. 



q,g 

f{g + q)^ -i[f7cos(f))-sin(f])] and f {g - q) 

When fi ~ 2mT, the wave function coefficients in 
Eq.dlTl reduce to 



depend on n as (sec notations in text) 



(82) 



T- 



2ne 



°° r^ j-P" sin [(e + e') nt\ 

T^-' ./I=oA=o4=o « + «' 

1 fp l\ 1 fp ~l 



p + l 



P-'I 



- 1 - 



where "— " and "+" are corresponding to even and odd 
integer of 0/27r, respectively. Numerical results show 
that Acr-^'' {n = 2mr) w -Acr'='' [n = (2n + 1) tt] w ^. 

When 51 — (2n ± 1/2) tt, wave function coefficients in 
Ea.(|47|) can be rewritten as 



E: 



e; 



(-l)"e-'^/4 _i 
p 2 ^P 

^(-l)"e-''^/4 



1 



(83) 



d; = -z^ 



2 ^'r> 

-^■Dl^-Zr, 
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and 



7 ('-n"p»'^/4 
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respectively. Hence one of products of matrix elements 
vanish V,;+'*'^j+-"'*-hV,-+'''^j+-"'' = 0. That renders the 
difference of conductivity independent on O ; 



Aa"" = t2- 



(=0 Jfe=0 Jp=0 



P^ sin [(£ + d) ^t\ ff __3 



(85) 



Et 



e; 

Dt 



2 

2 
-Dt = 1. 



{zp + 1) ; 
{zp ~ 1) 



Similarly for J7 = (2n -t- l)7r,the wave function coeffi- 
cients in Ea. (|T7)) become 



E' 



(-l)"(l-^) 



{Zp 



p p 

Thus the difference in conductivity due to barrier, 
Ea.([771). when fl is even or odd integer of tt, will clearly 
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